Comment on "Collective dynamics in liquid lithium, sodium, and aluminum" 
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In a recent paper, S. Singh and K. Tankeshwar (ST), [Phys. Rev. E 67, 012201 (2003)], proposed 
a new interpretation of the coUective dynamics in liquid metals, and, in particular, of the relaxation 
mechanisms ruling the density fluctuations propagation. At variance with both the predictions of 
the current literature and the results of recent Inelastic X-ray Scattering (IXS) experiments, ST 
associate the quasielastic component of the S{Q,lj) to the thermal relaxation, as it holds in an 
ordinary adiabatic hydrodynamics valid for non-conductive liquids and in the Q -^ limit. We 
show here that this interpretation leads to a non-physical behaviour of different thermodynamic 
and transport parameters. 

PACS numbers: 67.55.Jd; 67.40.Fd; 61.10.Eq; 63.50.-l-x 



Recently, Singh and Tankeshwar (ST) reported a new 
version of the finite wavevector, Q, extension of the gen- 
eraUzed hydrodynamics (the so called " molecular hydro- 
dynamics" ) applied to the dynamics of simple liquid met- 
als [l|. The goal of the ST paper is twofold. On one 
side, the authors propose an analytic expression for the 
density-density correlation function, F{Q,t), which is 
based on the existence of two relaxation channels, and 
they test it through its Fourier transform against the 
dynamic structure factor recently measured by inelastic 
X-ray scattering (IXS) on several liquid metals. On the 
other side, the authors assign a specific meaning to these 
two relaxation processes: the slower process is traced 
back to the thermal relaxation, while the faster one is 
ascribed to viscous effects. More specifically, at vari- 
ance with the results of previous studies performed on 
the same IXS data used by ST Q, y, Q , they assign the 
quasielastic component of the IXS spectra to the thermal 
relaxation. 

The aim of this comment is to point out how, although 
the analytical spectral shape proposed by the authors 
turns out to reproduce the experimental spectra at a 
reasonable level of accuracy, the interpretation behind 
it is inadequate to describe the microdynamics of sim- 
ple metals. In particular, assigning the whole quasielas- 
tic IXS spectrum to the thermal relaxation, the value 
and the Q-dependence of the derived fitting parameters 
are not consistent with independent determination of the 
transport (viscosity, thermal conductivity) and thermo- 
dynamic (specific heats) properties and disagrees with 
what it is expected on a general ground. 

The dynamic structure factor of simple liquids, in the 
Q —!■ limit, is constituted by three lines, a quasielastic, 
lorenzian, component centered at zero frequency, and a 
doublet (Brillouin peaks) symmetrically shifted at finite 
energy pj. In this limit, the energy position of the side 
lines, their broadening and the broadening of the cen- 
tral line are successfully described within the well known 



simple hydrodynamics approach. In particular, the fre- 
quency of the Brillouin mode is LUg — CqQ, being cq the 
adiabatic sound velocity, while its linewidth (HWHM) is 
ruled by thermal {Tth) and viscous (F,;) damping as 



= F, 



r*,, 



= C[(7-1)A/C, 



(1) 



where A, 7 and Cp are the thermal conductivity, specific 
heat ratio and constant pressure heat capacity, while 77^ 
and rjB are the shear and bulk viscosity, respectively. The 
width of the quasielastic mode, instead, is l/rth — DtQ^ , 
i.e. it is ruled by the thermal diffusion only. 

In Ref.yJl, ST extend this approach to the finite 
wavevector domain probed by IXS, and they apply it 
to a specific sub-class of simple liquids, namely liquid 
metals. To this purpose, ST replace the lorenzian shape 
stemming from simple hydrodynamics with a hyperbolic 
secant shape, leaving unchanged the origin of the indi- 
vidual contributions. 

The extension to finite wavevectors, however, in the 
special case of highly conductive systems, requires careful 
evaluation of the physics behind the model, for two main 
reasons. First, as pointed out in the celebrated textbook 
by Faber 'q, in a liquid metal, owing to the high thermal 
conductivity, on pushing the wavevector to values com- 
parable to the inverse mean interparticle distance (as in 
the present case) the quantity DtQ^ soon becomes larger 
than the Brillouin frequency ujb- Consequently, as pre- 
dicted (ref. |y|), the thermal peak broadens ultimately 
overlapping with the Brillouin lines, the sound propaga- 
tion turns from adiabatic to isothermal, and independent 
thermal fluctuations become impossible. A second rea- 
son is that, on the snapshot timescale probed by IXS 
(Thz), the diffusive atomic motion caracteristic of the 
liquid state looks frozen, and therefore part of the vis- 
cous contribution ruling the Brillouin width at low Q is 



transferred to the elastic line. Both of these reasons have 
been disregarded by ST. 

A formal way to quantify the scenario qualitatively de- 
picted above, is to introduce the memory function for- 
malism. Within this framework, the normalized density 
autocorrelation function 4'{Q,t) = F{Q,t)/S{Q) obeys 
the Langevin equation 



'(P{Q,t) + Lu^iQ)q^{Q,t)+ / M{Q,t-t')^{Q,t')dt' ^ 



where lo^{Q) = - ^ (Q,0) = KBTQ^/mS{Q) is related 
to the generalized isothermal sound velocity Ct{Q) = 
^o{Q)/Q and M{Q,t) is the second order memory func- 
tion which, in the hydrodynamic limit reads: 



M{Q, t) = 2DvQ^S{t) + (7 - 1) Wo(Q)e 



-DrQ^t 



(3) 



with Dy — {^Vs + Vb)/p the longitudinal viscosity and 7 
the specific heat ratio. It is easy to show from Eq. |j2Jl and 
©j indeed, that if the condition uig >> DxQ^ holds, 
the dynamic structure factor (i.e. the Fourier trasform 
of F) reduces to the sum of a (thermal) central loren- 
tian contribution and a damped harmonic oscillator spec- 
trum of characteristic frequency lob ~ -y/7'^o (which is 
the sum of two symmetrically shifted lorenzians in the 
limit r << ujq). As 7 = [ci^/ctf, one ultimately gets the 
hydrodynamic result, i.e. an adiabatic regime. 

The finite Q wavevector generalization stems dropping 
the hypothesis of the instantaneous (Markovian) nature 
of the viscous term. In particular, one can introduce a 
finite timescale (r) for the decay of the first term in Eq. 
(PI and, according to the ujqt value, the viscous relax- 
ation can affect both the quasielastic and the inelastic 
peaks. Actually, the situation is even more involved. In 
earlier MD studies of liquid metals it has been soon real- 
ized that the viscous dynamics in the microscopic regime 
(i.e. at wavelength comparable with the inverse mean 
inter-particle separation) proceeds through two distinct 
processes, characterized by two well separate time scales 
ija and T^) H- This idea has been recently substantiated 
by a number of IXS investigation, where the presence and 
the role of the two relaxation mechanisms has been iden- 
tified and widely discussed QiS Sialyl- In particular, it 
has been pointed out that the slower (-a) relaxation time 
satisfies the condition LOB{Q)Ta{Q) >> 1, i.e. some part 
of the viscous flow is frozen. As a consequence, at the 
wavevectors typical of the IXS experiments (Q = 1 ^ 20 
nm~^) the quasielastic spectrum acquires a component 
arising from this frozen structural relaxation. Moreover, 
owing to the high thermal conductivity, in liquid metals 
the condition LUBTth < 1 holds down to wavevectors of the 
order of 0.1 nm^^ 15]. In this limit, the sound propaga- 
tion is isothermal and not adiabatic, which means that 
the thermal relaxation is too broad to give a noticeable 



quasielastic component, while it renormalizes the sound 
velocity from the hydrodynamic value cq to the isother- 
mal value Cj, as clearly discussed in Ref. ^ and pointed 
out in previous literature jlfif . 

Summing up, the quasielastic peak in the S{Q,uj) is ex- 
pected to arise from the frozen structural relaxation and 
not from the thermal process. Its linewidth is therefore 
associated to t~^ and not to DtQ^- 

A decisive support to the origin that we propose for 
0(2)the quasielastic peak can be gained by the wavevector 
and temperature dependence of the spectrum. The sig- 
nature of a viscous origin of the central peak are, indeed, 
i) a substantial Q-independence (tq, is almost Q inde- 
pendent) , and ii) a strong temperature dependence of its 
width (the structural relaxation follows the Arrhenius or 
an even faster behaviour [l(J|). On the contrary, a central 
peak of thermal origin must show i) a strong Q depen- 
dence (width ex Q^ until the structure factor is flat) and 
ii) a very weak temperature dependence (traced back to 
the mild temperature dependence of the thermal conduc- 
tivity) . While the temperature dependence of the central 
peak is difficult to determine experimentally, because of 
the impossibility to supercool liquid metals (some hints 
could came from the molecular dynamics simulation, see 
below), the Q dependence of the central peak is rather 
clear: its experimental width (as determined by IXS spec- 
tra) is a quantity that can be roughly determined by a 
ruler, or better by a fitting procedure, and turns out to 
be almost constant in the examined Q-range. If one in- 
terpreted (following ST) this width as due to the thermal 
relaxation process, it would be equal to DtQ^ and, as a 
consequence, the generalized thermal diffusivity, would 
obviously become strongly Q-dependent ( Dt{Q) oc Q~^ 
) , ultimately dropping by a factor « 400 from the hydro- 
dynamic limit to the examined Q range. This behaviour 
is difficult to understand on a physical ground, and is very 
different from that observed in other liquids. Indeed, al- 
though we are not aware of any explicit Dt{Q) estimates 
for liquid metals, detailed calculations for other liquids 
(water jll|) show a decrease of only a factor ten on going 
from Q=0 and the Q value of the the maximum of the 
static structure factor. The too strong Q dependence of 
Dt{Q) derived by ST speaks against the assignment of 
the central peak to the thermal process. Turning now 
our attention to the temperature dependence of the cen- 
tral peak width, a recent MD work on underccoled liquid 
lithium showed that this width is actually strongly tem- 
perature dependent. More important, the derived relax- 
ation time. Ton closely follows the behaviour of the mass 
diffusion coefficient (as expected, for example, from the 
predictions of the Mode Coupling Theory for the struc- 
tural relaxation time |l2l|'). increasing by more than a 
factor ten in the spanned temperature range [iJl- Such 
a strong temperature dependence is not expected for the 
thermal diffusivity, indicating once more the non-thermal 
origin of the central peak. 
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FIG. 1: Comparison between the generalized specific heat 
ratio 7 = c-pjc^ as obtained by the ST model (o) and by 
direct numerical calculation (•) 



A further consequence of the assignment of the central 
peak to the thermal relaxation process is the anomalous 
Q dependence, obtained by ST, of the parameter describ- 
ing the ratio of the elastic to inelastic scattering intensity. 
In the authors notation this is the quantity a(ff) entering 
in their Eq. (4), which is related to the generalized spe- 
cific heat ratio 7(Q) via the expression 7 = j^^. In FigQ] 
we report 7(Q) as obtained by ST compared to previous 
numerical calculation for the case of liquid lithium 1^ . 
The values from the ST analysis clearly overestimates the 
calculated ones, although these latter have been obtained 
through simulations which well reproduce the experimen- 
tal spectra. 

Finally, the hydrodynamic expression (Eq. (Q) which 
ST utilize to estimate the sound attenuation (and to 
judge the reliability of their model) is valid only for or- 
dinary (non conducting) liquids in the Q — > limit. It 
accounts, indeed, for the first-order non vanishing atten- 
uation contribution from the thermal diffusive mode, and 
it relates the attenuation to the whole longitudinal vis- 
cosity. As pointed out in Ref. p], the small attenua- 
tion contribution from the thermal process in the IXS 
wavevectors range has, in liquid metals, the different ex- 
pression 



Vtn^(n-\)c\lDT 



(4) 



being q the isothermal sound speed. More important, at 
the probed wavevectors, only the microscopic, relaxed, 
part of the viscosity is involved in the sound damping 
F^, as the structural relaxation contribution if frozen on 
the probed time-scale (in liquid lithium this latter contri- 
bution from the microscopic part is one half of the whole 
viscous term |g|). 

In conclusion, in our opinion the model proposed by 



Singh and Tankeswar has the merit to indicate a possible 
heuristic analytical time- (and frequency-) dependence of 
the intermediate scattering function (dynamic structure 
factor) appropriate to catch some features of the dynam- 
ics in the Thz frequency region. Indeed, for instance, the 
functional form proposed by ST has all the spectral mo- 
ment finite, while, on the contrary, the memory function 
approach, leading to spectral shape that are described 
by fraction of polynoms, has only a finite number of con- 
verging spectral moments. The ST analysis, however, 
seems to be rather weak in the physical interpretation of 
the relaxation processes that drive the collective dynam- 
ics. None of the two parameters ti and T2 has a clear 
meaning and can be directly associated to any relaxation 
mechanism, as it happens for example within a memory 
function framework. As a consequence of the incorrect 
identification of the role of the thermal process, the pa- 
rameter T\ has no connection with any thermal property 
(being in fact mainly related to the structural relaxation 
time Tq), the generalized specific heat ratio obtained from 
the fitting shows marked discrepancies with the expected 
behavior, and the viscosity can not be estimated through 
Eq. (11) in Ref. [1]. 

We thank D. Fioretto, J. P. Hansen, F. Sacchetti and 
G. Senatore for stimulating discussions. Fruitful email 
exchange with S. Singh and K. Tankeshwar is also ac- 
knowledged. 
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